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The CMS Collaboration has released the results of its search for supersymmetry, by applying an 
aj method to 1.1/fb of data at 7TeV. The null result excludes (at 95% C.L.) a low-mass region 
of the Constrained MSSM's parameter space that was previously favored by other experiments. 
Additionally the negative result of the XENONioo dark matter search has excluded (at 90% C.L.) 
values of the spin-independent scattering cross sections o^^ as low as 10^^ pb. We incorporate 
these improved experimental constraints into a global Bayesian fit of the Constrained MSSM by 
constructing approximate likelihood functions. In the case of the ocj limit, we simulate detector 
efficiency for the CMS ccj 1.1/fb analysis and validate our method against the official 95% C.L. 
contour. We identify the 68% and 95% credible posterior regions of the CMSSM parameters, and 
also find the best-fit point. We find that the credible regions change considerably once a likelihood 
from a J- is included, in particular the narrow light Higgs resonance region becomes excluded, but 
the focus point/horizontal branch region remains allowed at the la level. Adding the limit from 
XENONioo has a weaker additional effect, in part due to large uncertainties in evaluating £7^', which 
we include in a conservative way, although we find that it reduces the posterior probability of the 
focus point region to the 2a level. The new regions of high posterior favor squarks lighter than the 
gluino and all but one Higgs bosons heavy. The dark matter neutralino mass is found in the range 
250 GeV < m^^ < 343 GeV (at la) while, as the result of improved limits from the LHC, the favored 

range of it^' is pushed down to values below 10^' pb. We highlight tension between S{g — 2)^,^^^ 
and BTZ (B — ^ Xgy), which is exacerbated by including the aj limit; each constraint favors a different 
region of the CMSSM's mass parameters. 
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I. INTRODUCTION 

The search for new physics at the Large Hadron Collider (LHC) began in earnest last year. An initial dataset of 
about 35 /pb at ^/s = 7TeV was followed this year by a much larger collection of about 5/fb of data. Based on 
about 1 /fb of analyzed data, earlier this year new limits were published by ATLAS fr] and CMS (2] experimental 
collaborations which significantly improved early LEP [3] and recent Tevatron results [4] as well as their own initial 
exclusion limits on the mass scale of low-energy supersymmetry (SUSY), in particular, on the mass parameter space 
of the Constrained MSSM (CMSSM) [5]. 

The CMSSM is a tractable unified model of effective softly broken supersymmetry, which includes models like, 
e.g., a relaxed version of the minimal supergravity model [6]. Despite its restrictive boundary conditions, the 
CMSSM has a rich phenomenology and has been extensively used as a basis for evaluating prospects for SUSY 
searches at the LHC and in other collider, noncollider and dark matter (DM) experiments. The model is defined by 
four continuous parameters and one sign [s]: niQ, the universal scalar mass; mi/2, the universal gaugino mass; Aq, 
the universal trilinear coupling; tan^, the ratio of the Higgs vacuum expectation values; and sgnji; the sign of the 
Higgs/Higgsino mass parameter p. We denote them collectively by 9. 

So far, the best limits from the LHC experiments ATLAS and CMS have come from channels involving only 
hadronic final states, with ATLAS investigating jets plus missing transverse momentum, and CMS using jets plus 
missing transverse energy, a "razor" and search methods. In particular, they involve different ways of 

clustering events into effective dijet systems. Based on 1.1 /fb of data the strongest bounds on the scale of SUSY 
have so far come from the ctj method. 

Even before the turn-on of the LHC, the parameter space of the CMSSM had been significantly constrained 0- 
[14!] by a variety of experimental data, most notably by LEP bounds on electroweak observables, masses of the 
Standard Model (SM) and SM like Higgs boson h and the lighter chargino l^]/ data on heavy flavour processes: 
S — )• Xs7, Bj, — )• TV and Bs — )• ji^n" , the difference between experimental and Standard Model contributions to 
the anomalous magnetic moment of the muon 5{g — 2)^^^^ [15], as well as the relic abundance Cl^h'^ Ii6il of the 
lightest neutralino, which is assumed to be the dominant component of cold dark matter in the Universe. 

In particular, a method favored a region of the CMSSM's parameter space with mi/2 ^ a few hundred GeV 
and niQ < ff2i/2 llial . In a Bayesian approach, this region was also favored but in addition a region of parameter 
space with substantial posterior probability was also found in the so-called horizontal branch, or focus point (1^, 
region of large niQ > 1 TeV and ^ niQ and in, partly overlapping with it, light Higgs resonance region [9, 14,]. 
A region of sizable posterior probability also exists in between those two "islands" of the highest probability. 

Whilst the first results from the LHC had a fairly mild effect on both regions fi^-20], the current limits based 
on the early 2011 dataset of 1/fb published so far take a much deeper bite into the {niQ, mi/2) plane of the 
CMSSM l2dEl- 

Searches for signals of DM in direct detection (DD) experiments have over the last few years also led to much 
improved limits [22]. Most notably, last Spring, XENONioo [23, 24J, with 100.9 days of data, excluded (at 90% C.L.) 
spin-independent (SI) scattering cross sections (7^^ as low as 10~^ pb, for some dark matter particle masses (23!] . The 
impact on the CMSSM's parameter space of this new XENONioo limit has been investigated in combination with 
2010 LHC limits in Refs. 1 20I, H, I26I1 and in combination with current LHC limits in Ref. [21.] . Moreover, LHCb has 
recently reported an upper bound on BTZ {Bg (13] that is smaller than the previous best upper bound. 

In this paper, our goal is to perform a Bayesian analysis of the CMSSM's parameter space that includes the 
currently strongest LHC limits on the SUSY mass scale and the XENONioo limit on the DM scattering cross 
section by carefully estimating associated uncertainties and including them in the likelihood function. Combining 
limits from different search channels at the LHC is a rather challenging task. In our approach we focus on the 
result derived by CMS from the a 7 method applied to 1.1/fb of data, since it currently provides the strongest 
constraint on the {iiiq, WI1/2) plane of the CMSSM (i^. EE El- Ell- will simulate the ocj experiment at the event 
level, assume that the experiment was a Poisson process and construct a likelihood map on the CMSSM's (otq, 
nil 12) plane. Our approach is similar to that of Ref. |iq]; however, it differs from that of Ref. fi^ I21I, I25II , in which 
the likelihood is modeled with an empirical formula, and from that of Ref. [20], in which the likelihood for the 
CMS limit is approximated with a step function. In the latter approaches, the likelihood is constructed from the 
published CMS olj 95% contour. In contrast, in our approach, the likelihood is constructed in the whole (mg, mi/2) 
plane and next validated against the official CMS a 7 95% contour. 

This paper is organized as follows. In SeclUwe detail our methodology, including our statistical tools, scanning 
algorithm, and our treatment of the likelihood from the CMS ctj 1.1/fb analysis. In Sec. |III]we present results of 
scans that include likelihoods from the CMS ocj 1.1/fb analysis and from XENONioo. We summarise our findings 
inSecHVl 
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II. METHOD 
A. The statistical framework 

Our goal is to identify the regions of the CMSSM's parameter space that are in best agreement with all relevant 
experimental constraints, including the aj and XENONioo limits. The mass spectra and other observables are also 

dependent on Standard Model parameters, most notably the top pole mass m^"^'^, the bottom mass mi,{mi,)^^ , the 

strong coupling as{Mz)^^ , and l/dgiaiMz)^^ ■ These "nuisance" parameters, which we collectively denote by cp, 
have been shown to play a significant role in a statistical treatment |8l4mli4'.'2q|. 

To set the stage, we define the best-fitting regions with Bayesian and frequentist statistics. In both approaches 
one considers the likelihood - the probability of obtaining experimental data for observables given the CMSSM's 
underlying parameters, 

c{e,cp)^p{d\e,cp). (1) 

To see the likelihood's dependence on a particular parameter or set of parameters, one maximizes over the CMSSM's 
other parameters, to obtain the profile likelihood, 

C{0) =rmixp{d\9,(p). (2) 

'P 

In Bayesian statistics one addresses the question of the posterior - what is the probability of the CMSSM's 
parameter values given the experimental data? To this end one employs Bayes's theorem to find the posterior 
probability density function (pdf); 

,fe.m - ^MMM. ,3) 

The Bayesian approach requires that we articulate our prior knowledge of the CMSSM's parameters in the prior, 
n{6,(p). Finally, the denominator, p{d) is the evidence, which, because we are not interested in model comparison, 
is merely a normalization factor. 

To see the posterior 's dependence on a particular parameter, or set of parameters, one integrates, or marginalizes, 
over the CMSSM's other parameters, as well as SM nuisance parameters, to obtain the marginalized posterior pdf. 

The two-dimensional region of the CMSSM's, e.g., parameter space (otq, ^1/2) that is in best agreement with the 
experiments, with respect to the posterior - the credible region - is the smallest region, R, that contains a given 
fraction of the posterior, that is, the smallest region such that 

Jp{mo,mi/2\d)dmodmi/2 = l-e. (4) 

R 

The one-dimensional credible region, however, is the region such that the posterior probability of the parameter 
being above the region is equal to the posterior probability of the parameter being below the region and equal to a 
given fraction of the posterior. That is, the one-dimensional credible region from L to LI of, e.g., wzq, satisfies 

L 00 

p{mo\d) dmo = J p(mo|d) dmo = -e. (5) 
u 

In the frequentist approach, the fc-dimensional region of the CMSSM's parameter space that is in best agreement 
with the experiments, with respect to the likelihood - the confidence interval - is the region in which the is 
within Ax^ of the minimum x^/ where Ax^ is such that 

F(Ax\k)=l-e, (6) 



and F (^Ax^, fc) is a cumulative x^ distribution with k degrees of freedom. 

The parameter point with the minimum x^ (or, equivalently, with the maximum likelihood) is the best-fit point. 
In frequentist statistics, the best-fit point has a special significance; the confidence intervals are constructed from 
the best-fit point. In contrast, the best-fit point has no significance in Bayesian statistics. 
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Parameter 


Description 


Prior Range 


Scale 


CMSSM 


mo 


Universal scalar mass 


100, 2000 


Log 


mi/2 


Universal gaugino mass 


100, 1000 


Log 


Ao 


Universal trilinear coupling 


-2000, 2000 


Linear 


tan^ 


Ratio of Higgs VEVs 


3, 62 


Linear 


sgnji 


Sign of Higgs parameter 


+1 


Fixed 


Nuisance 


pole 


Top quark pole mass 


163.7, 178.1 


Linear 




Bottom quark mass 


3.92, 4.48 


Linear 




Strong coupling 


0.1096, 0.1256 


Linear 


l/aem(Mz)'^S 


Reciprocal of electromagnetic coupling 


127.846, 127.99 


Linear 



Table I: Priors for the CMSSM's parameters and for the Standard Model's nuisance parameters that we used in our 

scans. Masses are in GeV. 



From the best-fit point, one can obtain a p-value: the probability of obtaining a value from experimental 
measurements equal or larger than the best-fit x^/ accounting for the number of degrees of freedom in the fit, 

p-value = 1 - f (^x^, nj , (7) 

where n, the number of degrees of freedom in the fit, is the number of experimental constraints in the calculation 
minus the number of model parameters that were fitted. 

Our credible regions and confidence intervals are defined by Eq. |4| Eq. [5] and Eq. [6l with e = 0.32 for la and 
e = 0.05 for la. In Eq. [6l for a two-dimensional confidence interval, this corresponds to Ax^ = 2.30 for Icr and 
Ax^ = 5.99 for 2a. 

Scanning the CMSSM's parameter space is computationally intensive - a simple grid-scan of an eight-dimensional 
space is impractical. To scan the CMSSM's parameter space efficiently, we use a modern (2006) Monte Carlo al- 
gorithm, called Nested Sampling [^0], which is tailored to work with Bayesian statistics. We chose the Nested 
Sampling settings so that the algorithm would accurately map the posterior, but not necessarily the likelihood, in 
a reasonable CPU time (4000 live points and a stopping condition of 0.5). 

Because of our lack of prior knowledge of the CMSSM's parameters, we invoke the principle of insufficient 
reason and, as previously in Jot bji n^fl we choose noninformative priors for the CMSSM's parameters that 
equally weight either linear or logarithmic intervals. The priors that we choose for niQ and WI1/2 equally weight 
logarithmic intervals (log priors). This choice has been shown [ij not to suffer from the volume effect, unlike 
the flat prior, and it also reduces the amount of fine tuning needed to achieve radiative electroweak symmetry 
breaking. For tan j6, Aq and for the SM's nuisance parameters we choose equally weighted linear intervals (linear 
priors). The prior ranges of the CMSSM's parameters and of the SM's nuisance parameters over which we scan are 
listed in Table ID 

We use our updated and modified version of the SuperBayeS |@>E3] computer program to perform four scans of 
the CMSSM: 

1. To set the stage for examining the impact of LHC and XENONioo limits, a scan involving constraints from 
only non-LHC experiments, including those listed above (see Table HU] for a complete list of observables and 
their values) but without a likelihood from XENONioo. 

2. To validate our method of including the LHC constraints, a scan with a likelihood that we derived from the 
CMS a J 1.1 /fb analysis, and a likelihood from the experiments that constrain the Standard Model's nuisance 
parameters only. 

3. To examine the impact of the current LHC constraints, a scan with a likelihood from non-LHC experiments 
and from the ctj but without a likelihood from XENONioo. 

4. To see the additional impact of the current XENONioo limit scan with a likelihood from non-LHC experi- 
ments, the Ut and from XENONioo. 
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Ht Bin (GeV) 


275 - 325 


325 - 375 


375 - 475 


475 - 575 


575 - 675 


675 - 775 


775 - 875 


> 875 


pleadmg ^^^^^ 


73 


87 


100 


100 


100 


100 


100 


100 


psecond ((3gv) 


73 


87 


100 


100 


100 


100 


100 


100 


p?,*'^''^ (GeV) 


37 


43 


50 


50 


50 


50 


50 


50 


Observed events 


> 0.55 


782 


321 


196 


62 


21 


6 


3 


1 



Table II: Definition of the Hj bins and the corresponding pj thresholds for the leading, second, and all other 
remaining jets in the event. Observed events refers to the number of events passing all (Xj cuts for 1.1 /fb of data 

collected by the CMS Collaboration [3^. 



B. The efficiency maps for the CMS kj 1.1/fb analysis 

We derived our LHC likelihood for the CMS search (33'] for R-parity conserving supersymmetry in all-hadronic 
events via a kinematic variable Oij. The results based on the LHC data sample of 1.1/fb of integrated luminosity 
recorded at = 7 TeV showed no excess of events over the SM predictions. Our aim was to translate the analysis 
scheme into a simplified approach to obtain the signal selection efficiency for a large number of points in the 
CMSSM parameter space. To this end we generated a map of points for the CMSSM parameters mg in the range of 
(50, 2000) GeV and mj /2 in the range of (50, 1000) GeV, both with a step of 50 GeV. We fixed the values of the other 
CMSSM parameters: Aq = 0, tan /3 = 10 and sgn ]i = 1. For each point, for a defined set of the CMSSM parameters 
we calculated a mass spectrum and a table of supersymmetric particle decays using the programs SoftSUSY Q 
and SUSY-HIT [^1^], respectively. For each point, we then generated 10,000 events with the Monte Carlo generator 
PYTHIA [36] with a cross section value obtained at the leading order. We analyzed events at the generator level. 
Only the geometrical acceptance was applied to simulate the CMS detector response. Leptons were accepted in 
the pseudorapidity range rj < 2.5 for electrons and rj < 2.1 for muons, respectively. The isolation of leptons was 
checked at the generator level. Stable, generator level particles, excluding neutrinos, were clustered into jets using 
the anti-fcj algorithm with a cone size parameter R — 0.5. Jets were accepted up to |;/| < 3. 

The aim of the aj analysis was to select hadronic events with high transverse momentum jets. Initially, events 
with at least one jet with transverse momentum > 50 GeV and i] > 3 were accepted if no isolated lepton or 
photon were found in the event, namely events with an isolated lepton (electron or muon) with pj > 10 GeV or an 
isolated photon with pj > 25 GeV were rejected. Events had to satisfy a condition based on an Hj variable, defined 

as Hj = YJi=\ Ej^'; Hj was required to be above 275 GeV. Following the CMS analysis f^^, the trigger was fully 
efficient for selected events, and therefore no attempt to emulate the trigger was made. 

The offline analysis used the following Hj binning: 275 - 325, 325 - 375, 375 - 475, 475 - 575, 575 - 675, 675 - 
775' 775 " 875, and > 875 GeV. In each Hj bin, we applied the same cuts on transverse momentum of jets (leading, 
second, others) in the event as in Ref. [^^]. Details are shown in Table Ull 

The main discriminator against QCD multijet production is the variable a j defined for dijet events as 

^T=^= ^ (8) 

v/(ELEf)^-(ELpr')2-(ELpf)^ 

where E^j^ is the transverse energy of the less energetic jet in the event with two jets and Mj is a transverse mass 
of the dijet system defined above. 

In events with more than two jets, two pseudojets were formed following the same strategy as in Ref. in such 
a way that the Ej difference between two pseudojets was minimised. The value of £7 of each of the two pseudojets 

was obtained by a scalar summing of the contributing Ej* of jets. In the ideal case, the dijet system had Ej*^ = Ej*^ 
and jets are back-to-back which resulted in a limit of dlj = 0.5, where the momentum of jets is large compared to 

the masses of jets. For back-to-back jets with Ej^^ 7^ Ej'^, values of ocj were smaller than 0.5. Signal events with 
missing transverse energy resulted in xj greater than 0.5. Therefore the QCD background was efficiently rejected 
with the final cut of ccj > 0.55. The remaining events after all cuts were compared with SM background expectation 
predictions. No excess was found and the numbers listed in Table Ull were used as numbers of observed background 
events. 

For the selection defined above, we prepared eight efficiency maps, one for each Hj bin, two of which are shown 
as examples in Fig. [1] 

Experimental selections involving missing transverse energy estimates and selections from an analysis of the 
distance between a jet and Hj were not implemented, because the relevant variables calculated at the generator 
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(b) The case of 325 GeV < Hj < 375 GeV. 



Figure 1: The efficiency maps of our approximation to the CMS olj 1.1/fb analysis on the (otq, ^1/2) plane, for 

tan [i = 10 and Aq = 0, for the two dominant Hj bins. 



level are not reliable. Nevertheless, as discussed below, good agreement between efficiency maps obtained in this 
analysis, and the CMS result shows that such an approximation is justified. 



C. The likelihood from the CMS arl-l/fh analysis 

The approximate efficiency maps derived in Sec. lIIBI for the CMS ocj l-l/fb limit allow us to evaluate a likelihood, 
so that we can find the regions of the CMSSM's parameter space in best agreement with the currently strongest 
LHC limit. 

We assume that the experiment can be described by a Poisson distribution; that is, that the events were indepen- 
dent and that the likelihood of observing events was described by a Poisson distribution. We assume that the total 
number of observed events was the sum of contributions from supersymmetric processes and from SM processes. 
Following the dj cuts described in Sec. IIIBi we consider separately the eight bins of the kinematical variable Hj 
summarized in Table HH 

The number of supersymmetric events that we expected in an z'-th Hj bin, s,, is the product of the detector 
efficiency for that bin (the fraction of events that survive the ocj cuts), e,, the integrated luminosity, J L = 1.1/fb, 
and the total cross section for the production of supersymmetric particles at y/s — 7 TeV, a, 

Si = e,- X (7 X J L. (9) 

The likelihood for this case, £, which is the probability of observing a set of {0,} events given that we expected 
{s/} supersymmetric events and {fo;} SM background events, is a product of Poisson distribution for each bin with 
means A; = s,- + fo,-, 

c=n- — ^i±^. (10) 

We assume that Standard Model backgrounds, b,, are precisely known. If their experimental errors were significant, 
the expression for the likelihood would be multiplied by distributions describing the SM backgrounds [^8] and they 
would be included as nuisance parameters. This is not the case for the ocj analysis. 

All-hadronic final-state processes are expected to be independent of tan [i and Aq, because these parameters have 
little influence on the squark and gluino masses, and this has been confirmed by the CMS Collaboration in the 
CMS Xj 1.1/fb analysis, as well as, e.g., in Ref. [iQ, ^q], and also in this study. For this reason we created efficiency 
maps using the mentioned earlier fixed-grid scan with fixed values of tan^ = 10 and Aq = 0, as described 
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using Nested Sampling and a log prior 



Figure 2: The maps of (a) the likelihood and (b) the profile likelihood for our approximation to the CMS ccj 1.1/fb 
analysis on the {mq, mi/2) plane. The dashed blue lines show the official CMS 95% exclusion contour. 



in Sec. IIIBI The cross section was calculated to leading order with PYTHIA [^6]0 We used this information to 
produce a likelihood map using Eq. [9] and Eq. [lol 

The resulting likelihood map is shown in Fig. [211 It exhibits the correct behavior; below the official 95% contour 
obtained by CMS ocj 1.1/fb analysis 1I33I1 , which is also indicated, the likelihood is small. Our approximate 95% 
contour was calculated with the help of Eq. IH It corresponds to the boundary of the la range of the profile 
likelihood obtained using Nested Sampling and a log prior, which is shown in Fig. bbl and is very close to the 
official CMS contour which is also shown in the Figure. For each CMSSM point for which our scanning algorithm 
wished to evaluate the likelihood, we interpolated the likelihood from our likelihood map. 



III. RESULTS 



In this section we will present our numerical results. To start with, in Table H] we show prior ranges and dis- 
tributions of CMSSM parameters and of SM nuisance parameters, while in Table |III] we list the observables that 
we will use in our analysis. These come in three sets. First, by "Non-LHC" we collectively denote all relevant 
constraints from dark matter, O^h^, precision measurements: sin^flgff, M^j, etc, flavour physics: BTZ (B Xsj), 
BTZ {Bu ^ Tv), AMb,, and BTZ {Bs f(+F")EI the excess in the anomalous magnetic moment of the muon, S{g — 
2)^1^^^/ as well as LEP and Tevatron limits on the Higgs sector and superpartner masses. Second, we will apply cur- 
rently the best limit published by the CMS Collaboration, which is from its oij analysis of 1.1/fb of data. Finally, at 
the end we will also apply an upper limit on the elastic scattering of neutralino dark matter recently published 
by XENONioo, in order to investigate its additional impact on the CMSSM's parameters and various observables. 

Below we will present the results of our scans as one-dimensional (iD) or two-dimensional (2D) marginalized 
posterior pdf maps of the CMSSM's parameters and observables. In evaluating the posterior pdfs, we marginalized 
over the CMSSM's other parameters and the SM's nuisance parameters. 



' The cross section, and consequently the number of expected supersymmetric events, changes by over 10 orders of magnitude over the 
(mo,mi/2) plane. The resulting likelihood function is not, therefore, sensitive to next-to-leading order corrections to the cross section. Even if 
''NLO ~ fLO/ the corrections would only slightly shift the isocontours of the cross section and likelihood on the (mg, 'Wi/2) plane. 

^ Actually, LHCb I27], an LHC experiment, recently obtained the best upper limit of BTZ (Bs yi^ji^) < 1.5 x 10^*, but in this study we will 
nevertheless include it in the "Non-LHC" group of constraints. 
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Measurement 


Mean 


Exp. error 


The. error 1 Likelihood distribution 


Reference 


CMS aj- 1.1/fb analysis 




See text 


See text 





Poisson 




XENONioo 




< / {nij^), see text 





1000% 


Upper limit: error function 




Non-LHC 




0.1120 


0.0056 


10% 


Gaussian 


[1^ 


sin^ flgff 


0.231 16 


0.000 13 


0.000 15 


Gaussian 




Mw 


80.399 


0.023 


0.015 


Gaussian 


Uil 




30.5 


8.6 


1.0 


Gaussian 


Uil 


BTl (B ^ Xs7) X 10* 


3.60 


0.23 


0.21 


Gaussian 




BTZ (B„ ^ Tv) X 10^ 


1.66 


0.66 


0.38 


Gaussian 




4Mb, 


17.77 


0.12 


2.40 


Gaussian 






< 1.5 X 10-8 





14% 


Upper limit: error function 


[2Z] 


Nuisance 


l/aem(Mz)'^S 


127.916 


0.015 





Gaussian 






172.9 


1.1 





Gaussian 






4.19 


0.12 





Gaussian 






0.1184 


0.0006 





Gaussian 




LEP and Tevatron 95% limits 




> 114.4 





3 


Lower limit: error function 




a 


</K) 

> 50 








5% 


Upper limit - step function 
Lower limit: error function 


[42] ( \A^]) 


m ± 

Xi 

me-R 


> 103.5 (92.4) 





5% 


Lower limit: error function 


[^[ ( [4^) 


> 100 (73) 





5% 


Lower limit: error function 


[al([^) 




> 95 (73) 





5% 


Lower limit: error function 


[al([n]) 




> 87 (73) 





5% 


Lower limit: error function 


[al([^) 


m 


> 94 (43) 





5% 


Lower limit: error function 


[44] ( [1^]) 




> 95 (65) 





5% 


Lower limit: error function 


[^] ( [4^]) 




> 95 (59) 





5% 


Lower limit: error function 






> 375 





5% 


Lower limit: error function 


[4] 


nig 


> 289 





5% 


Lower limit: error function 


[4] 



Table III: The experimental measurements that constrain the CMSSM's parameters and the Standard Model's 
nuisance parameters. Masses are in GeV. The numbers in parentheses in the list of LEP and Tevatron experimental 
measurements are weaker experimental bounds, which we use for some sparticle mass hierarchies. 



A. Impact of the « j limit 

First, in Fig. [3} on the (mg, nii/2) plane, we show the results of a scan with a likelihood from the non-LHC 
constraints (left panel) and with the additional impact of imposing the a 7 1.1 /fb constraint (right panel). Dark 
(light) blue regions denote the la {2a) posterior pdf regions of the {wiq, mi/2) plane. In Fig. [33] one can see two 
distinct modes on the {niQ, mi/2) plane. The vertical la mode is located in in the stau coannihilation/ A-funnel 
(SC/AF) region, while the associated 2a region is caused by the AF only. On the opposite side of the diagonal of 
the {niQ, mi/2) plane, the horizontal mode corresponds to overlapping contributions from the light Higgs resonance 
region (the narrow la strip) and, above it, the 2a focus point (FP) /horizontal branch (HB) region. 

As a result, the posterior mean (denoted by a solid black dot) lies between the two modes, but closer to the SC/AF 
region, where the best-fit point (denoted by an encircled cross) is also located. The blue dashed line denoting the 
95% lower limit derived by the CMS ocj 1.1/ analysis is marked, but not applied here. 

The existence of the two broad regions is consistent with the findings of previous pre-LHC Bayesian analyzes |Q, 
[iil, [i4l[3ill , although their relative size does depend on the choice of the prior. Physically they both arise as a result 
of reducing the relic density of the neutralino, which is generally too large in the part of the (ffig, nii 12) plane where 
X is the lightest superpartner. In the SC region the density is reduced by a coannihilation with the lighter stau (and 
other sleptons), and likewise in the AF region it is reduced by neutralino pair annihilation through the pseudoscalar 
Higgs A-funnel, hence the name. Likewise, the same mechanism is at play in the narrow horizontal light Higgs 
resonance region. In the FP region, on the other hand, as one moves down along a line roughly perpendicular 
to the diagonal of the (mo/ '^1/2) plane, the value of drops down from large values (for which the Lightest 
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Figure 3: Marginalised posterior pdf on the (mq, mi/2) plane, [(a)| before and |(b)| after we included a likelihood 

from the CMS a 7 1.1/fb limit. 



Supersjrmmetric Particle (LSP) is very binolike, with too large O^h^), and eventually becomes negative, implying 
a failure of radiative electroweak symmetry breaking conditions. Close to this boundary, in a rather narrow strip 
of the plane, }i is of order nii /2, the LSP develops a sizable Higgsino component, and the relic abundance becomes 
acceptable. 

The impact of the CMS ocj 1.1/fb limit, which is implemented in our analysis by simulating detector efficiency 
and evaluating the corresponding likelihood, as described in Sec. IlICi is shown in Fig. [3b] Clearly, the ocj limit 
has cut deep into the (mg, mi/2) plane's high posterior probability regions - a nominal fraction of the posterior pdf 
lu region is outside of the ctj 95% confidence interval. The ecj limit pushed the credible regions, as well as the 
best-fit point, to larger values of ?Mi/2- Significantly, the two modes on the {niQ, mi/2) plane have remained. The la 
SC / AF region has been pushed up nearly vertically, while the 2a one of AF only has become inflated and extended 
to larger values of (mg, mi/2)- Note also that the horizontal light Higgs resonance region has now completely 
disappeared. On the other hand, the FP/HB region is not excluded by the CMS aj 1.1 / fb limit; rather, it is pushed 
to larger values of niQ and, to a lesser extent, mi/2- 

An analogous comparison of the pre- and post-aj situation in the (Ag, tanj6) plane is presented in Fig. [4a] and 
Fig. 4b There are again two lu modes: one at relatively small values of tan /5 < 20 and the other one at much larger 
values 55), although at 2a an entire range of scanned values of tan /5 is allowed. The first mode corresponds 
to the la stau coannihilation region in Fig. |3b| however, in the much larger, 2a, mostly A-funnel region of the 
(mg,mi/2) plane much larger values of tan fi ^ 55 are predominant. On the other hand, in in the focus point region 
of large mg in Fig. [3b] we find 25 < tan j6 < 55. 

We can see that the application of the ccj constraint narrows down both modes in the posterior pdf on the 
(Ag, tan /5) plane, primarily as a result of pushing up and out the focus point region, with the effect of strongly 
disfavoring midrange values of tan/5. On the other hand the CMS limit shows fairly little effect on Ag, which 
remains poorly determined. 

Although our aj likelihood is independent of tan /5 and Ag, it is clear from Fig. [4] that ctj does impact these 
parameters. The CMSSM's parameters are correlated in a nontrivial way; once the ecj limit is added to the likelihood 
it alters the best-fitting values of mg and mi/2 and the values of tan /3 and Ag must be adjusted accordingly, so that 
the predicted values for the CMSSM's observables maintain agreement with the experimental measurements. 

In Fig.|5]we show iD marginalized posterior pdf plots for the CMSSM's parameters, constrained by the non-LHC 
experiments and by the a 7- limit. We also show the dark (light) blue horizontal bars which indicate the one- (two-) 
dimensional central credible regions. The high probability modes in the mg distributions at mg ^ 200 GeV in Fig. [3a] 
and in the distribution at mi/2 ^ 700 GeV in Fig. |5b| correspond to the high probability mode in the SC/AF 
region in Fig.|3b] Fig. |5c] shows that the experimental constraints favor positive values of Ag. The two modes in the 
tan /5 distribution in Fig. |3d| correspond to the two modes in Fig. [4] 

After marginalisation, it is difficult to identify in the mg iD marginalized posterior pdf the FP/HB mode, which 



10 



BayesFiTS (2011) 



BayesFiTS {2011) 



CO. 
C 



60 
50 
40 
30 



20 - 



10 - 



Posterior pdf 
CMSSM, ;(> 
Log priors 
Non-LHC 



lltr region 
2(7 region 
► Posterior mean 
} Best fit 



f 




-2-1012 

Aq (TeV) 

(a) The CMSSM's parameters constrained by the 
non-LHC experiments. 



to 



60 
50 
40 
30 
20 
10 



Posterior pdf 


' . 1 ' 
111 1(7 region 


■ CMSSM, }i>0 


2e region 


. Log priors 


• TOsterior meaxi 


Non-LHC + ffT 1.1 fb" 


' ® Best fit 




• 




9 



-2 



-1 



12 

Ao (TeV) 

(b) The CMSSM's parameters constrained by the 
CMS cij 1.1/fb analysis and non-LHC experiments. 



Figure 4: Marginalised posterior pdf on the (Aq, tan j5) plane, (a) before and (b) after we included a likelihood 

from the CMS a r 1.1/fb limit. 



was present in the 2D marginalized posterior pdf in Fig. [3} though it is present as a second smaller mode in the 
nil /2 iD marginalized posterior pdf in Fig. This is caused by the contribution of the large 2cr credible region at 
both large mi/2 and niQ. 

In Fig. [6] we present iD plots of marginalized posterior pdf for notable particle masses, when we scanned the 
CMSSM with a likelihood from aj and from the non-LHC experiments. These distributions are typically bimodal; 
one mode is from the mode in the SC/AF region and one mode is from the mode in the FP/HB region on the (mo, 
nii/i) plane in Fig. |3] 

Fig. |6a] shows that the posterior favors a moderately heavy lightest squarlH (the lighter stop) with ~ 500 GeV, 
corresponding to the mode in the SC/AF region with small mo ^ 100 GeV. There is a second, much smaller, mode 
in the posterior at niq ~ 1250 GeV, corresponding to the mode in the FP/HB region with large niQ ^ 2000 GeV. On 
the other hand, the squarks of the first two generations can be much heavier, up to ^ 2.7 TeV at la. 

Heavy gluinos, much heavier than lightest squarks, are favored with itig ~ 1500 GeV, as illustrated in Fig. |6bl 
This dominant range corresponds to the mode in the SC/AF region with large mi/2 ^ 700 GeV. Lighter gluinos 
with nig < 1000 GeV are not, however, excluded, by the credible region. This is a result of the mode in the FP/HB 
region. 

Fig. |6c] shows that the posterior pdf favors a lightest neutralino mass of ~ 300 GeV. This corresponds to a 
binoiike neutralino, with ^ QAnii/2- 

Lastly, Fig. I6dl shows the iD posterior pdf for the mass of the lightest Higgs in the CMSSM, in agreement with 
pre-LHC results 

The lightest Higgs boson in the CMSSM is to a very good approximation SM like and the LEP limit (114.4 GeV) 
applies. Note, however, that predicted Higgs masses below the LEP limit are permitted, because our likelihood 
function includes a 3 GeV theoretical error in the predicted Higgs mass. On the other side, the iD posterior pdf 
is rather narrow, with the 95% credible region in the range from 112.2 GeV to 119.2 GeV and the best-fit value of 
114.4 GeV. 

Being SM like, the lightest Higgs boson may well be the only Higgs state of the CMSSM accessible in Run 1. The 
other Higgses are typically nearly mass degenerate, and much heavier than the lightest Higgs. This is illustrated 
in Fig. [7| where we present a 2D posterior pdf map on the (m^, tan j6) plane in two cases: pre-a j (left panel) and 
post-aj (right panel). Fig. [t] shows that the aj limit favors a heavier CP-odd neutral Higgs, because it favors higher 
values of wlq. Intermediate values of tan^ ~ 30, however, are disfavored by a.j, and, consequently, very heavy 
CP-odd neutral Higgs with > 1500 GeV become disfavored. 



3 The mass of the lightest squark m,j = min {niq.) 
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Figure 5: One-dimensional marginalized posterior pdf for the CMSSM's parameters constrained by the 
CMS 1.1 /fb analysis and by the non-LHC experiments. The dark (light) blue horizontal bars span the one- 

(two-)dimensional central credible regions. 



B. Impact of the XENONioo limit 

The XENONioo Collaboration has recently published a new 90% C.L. exclusion limit on the {m^, cr^^) plane 
which significantly improves their previous preliminary limit {2^ . In this section we examine the impact of this 
new limit on the CMSSM's parameters and observables. 

A proper implementation of the XENONioo exclusion limit in our likelihood function is somewhat tricky. The 
90% exclusion contour alone is insufficient to reconstruct the likelihood function. On top of it, there are significant 
errors, which originate from poorly known inputs which are needed in evaluating a limit on from experimental 
data. First, one assumes a "default" value for the local density of 0.3 GeV/ c^, but errors in astrophysical parameters 
which are required to determine the local halo and its velocity distribution can result in errors in cr^^ of order 2 
although a much better determination has been claimed (IS]. More significantly, the fractional error in evaluating 
CTp^ coming from uncertainties in inputs to hadronic matrix elements can be of order 5, as discussed in detail in [2=;]. 
The main uncertainty is due to a poor knowledge of the 7T-nucleon a term, I^rcN/ which determines the strange 
quark component of the nucleon. 

Given these large uncertainties, we neglect the experimental error in the XENONioo 90% exclusion contour and 
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Figure 6: The mass of the lightest [(a)] squark, |(b)| gluino, (c) neutrahno, and (d) Higgs-boson in the CMSSM 
constrained by dj and by the non-LHC experiments. The horizontal bars have been defined in Fig. [5] 



approximate it in the likelihood with a step function, jC,{cr^^, m^;^-) = 0(1) for cr^^ > {<)(7'^^gQo/^{m^), where crp^gQo/^{m^) 
is the DM mass dependent XENONioo 90% limit. To incorporate the above errors in evaluating a^^ in a conservative 
way, we convolute the step function with a Gaussian with p = cr^^ and cr = 10 x a^^, resulting in a Gaussian error 
function. This will cause a rather large smearing out of the XENONioo limit. 

This can be seen in Fig. |8] where we present the impact of the XENONioo 90% C.L. limit (denoted with solid 
red curve) on probability maps of the {m^, (7^^) plane for our scans. In Fig. |8a]the XENONioo limit is not added 
to the likelihood, while in Fig. |8b] it is applied. By comparing both panels we can see that, once LHC limits have 
been applied (left panel), the la posterior region on the (;«;^, Up^) plane is only weakly affected by the additional 
XENONioo limit (right panel). 

Some interesting effects can nevertheless be noticed. First, a small 2a region above the XENONioo exclusion 
curve has shrunk somewhat, especially on the side of larger a^^, but, because of the large theoretical error assumed 
in this analysis, it remains allowed. Second, the large smearing affects the favored regions of also below the 
experimental curve. Note that, before the XENONioo is applied (Fig. l8al l there are actually two Icr regions close 
to each other. The lower one comes entirely from the stau-coannhilation region of small niQ and large fKi/2- The 
other one, just above it (along with a broader 2a region, both decreasing with m^), corresponds to the broad 2a AF 
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region in the (niQ, nii j^) plane. Once the XENONioo limit is added to the likelihood (Fig. ISbt , the largest values 
of (J^, just below the experimental curve, become excluded and the statistical significance of the whole AF region 
becomes reduced to the Icr level. On the other hand, the lower \(j region remains nearly intact. 

Clearly, recent LHC limits have had the effect of pushing down the most favored ranges of cr^, for the most part 
below ^ 10^^ pb, while pre-LHC data favored largest \a posterior ranges of cr^ at least an order of magnitude 
higher; compare, e.g., [fij. I31I] . This implies much poorer prospects for XENONioo, with expected reach of ^ 
10^^ pb to explore the SI cross sections currently favored by the CMSSM. 

It is also clear that, in light of the impact of LHC limits, one-ton detectors with planned sensitivity reach of 
~ 10"^*^ pb will now be needed to explore the most probable ranges of the (wi;^, <J^^ plane in the CMSSM. We also 
note that the most probable range of dark matter particle mass, 250 GeV < < 343 GeV (at Xcr) is somewhat above 
the range of the highest sensitivity of most detectors. Finally, in Table HVl we show the combined impact of the olj 
and xenon limits on the posterior ranges of several particle masses already constrained by other (Non-LHC) data. 
The Icr and Icr posterior ranges were calculated with Eq. [5] 

Finally, in Fig. [gj we show the impact of the XENONioo limit on the shapes of the 2D marginalized posterior 
pdf maps for the CMSSM parameters. The results shown here correspond to a scan with a likelihood from the 
non-LHC constraints, from the CMS dj 1.1/fb analysis and from XENONioo. The panels should be compared 
with the corresponding panels in Fig. [3] and Fig. [4] 

The discussion of Fig.|8]above helps in understanding new effects in Fig. [9] First, it is clear that, the configurations 
of the CMSSM parameters are already constrained by current LHC limits to give for the most part below the 
XENONioo limit. Unsurprisingly, adding it to the likelihood has a relatively weak additional impact on the 
CMSSM's parameters. The Xu posterior region in the high probability SC region on the (mg, m\ /2) plane in Fig. [ga] 
remains unaffected because the corresponding values of cr^ are the lowest. The change in the best-fit point is also 
probably not statistically significant. However, note that the broader 2a region of AF has shrunk somewhat on the 
side of large mg, especially in the direction of the HB/FP region. This is because in that direction the Higgsino 
component of the DM neutralino increases, causing in turn to increase and become at some point constrained 
by the XENONioo limit, especially for large tan j6. 

Moving farther toward the lower-right corner, we can see that the XENONioo result disfavors also the HB/FP 
region - the la credible region in the (mg, nii /2} plane in Fig. |3a] has now become only a 2a region. This is not 
surprising since predicted values of a^ in that regions are among the largest predicted in the CMSSM, again 
because of the increased Higgsino component of the neutralino. We stress, however, that the HB/FP region is not 
excluded by the XENONioo limit; rather, it is inside the 2a credible region. This follows from our conservative 
estimation of the error in the a'^ calculation. (We weakened the XENONioo limit by smearing it with a Gaussian 
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describing the theoretical error in the cr^ calculation.) We have checked that, if we assumed the theoretical error 
in cr^^ to be 0.1 x cr^, the XENONioo limit would exclude the focus point region at 2cr. Finally we note that the 
favored regions on the (Aq, tan j6) plane in Fig. 9b are not strongly affected by XENONioo, though the two modes 
on the (Ag, tan /5) plane are no longer connected. 



C. Prior dependence and the best-fit point 

We will now comment on the prior dependence of the results presented here, and will also come back to the 
discussion of the best-fit point. 

As is well known, some of the most challenging aspects of Bayesian statistics are the necessity to choose a prior 
and the sensitivity of the posterior to that choice. The prior dependence of the posterior is a measure of the 
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lack of the constraining power of the likelihood function. If the information from data included in the likelihood 
is sufficient to select ranges of model's parameters giving good fit to the constraints (high posterior probability 
regions), then the sensitivity of such ranges to the choice of priors should be weak, or even marginal. In such cases 
the likelihood has the ability to overpower noninformative priors, such as log priors and linear (flat) priors, and 
the resulting posterior has little prior dependence. If, however, the posterior is dependent on the choice of prior, 
the likelihood ought to be regarded as too weak to support robust conclusions about the model. 

In several pre-LHC studies prior dependence of the CMSSM was found to be substantial (see, e.g., fij,'^), and 
the situation in less constrained models was found to be even less satisfactory [47]. This was merely a reflection of 
poor constraining power of the data available at that time. It was also concluded that the choice of the prior that is 
linear in the log of the masses is more motivated by both physical and statistical reasons [14!,'^. From the physical 
point of view, log priors explore in much greater detail the low-mass region, where one typically needs less fine 
tuning in order to achieve radiative electroweak symmetry breaking. From the statistical point of view, log priors 
give the same a priori weight to all orders of magnitude in the masses, and thus appear to be less biased to giving 
larger statistical a priori weights to the large mass region, which under a flat prior has a much larger volume in 
parameter space. 

Motivated by the above arguments, in this study we have chosen log priors for the CMSSM's mass parameters, 
and our posterior distributions are, of course, dependent on the choice. As a way of examining the degree of the 
dependence, we identically repeated our scan of the CMSSM with a likelihood from the non-LHC experiments 
and from ecj, except that this time we chose linear priors for all the CMSSM's parameters. The resulting posterior 
distributions, which are shown in Fig. (To) are actually broadly similar to the distributions that resulted from log 
priors. With linear priors, we found three la modes in the marginalized posterior pdf on the (;«o, '^1/2) plane. 



rather than two as in Fig. 3b The third mode is present in Fig. [To] both large wJi/2 and wiq, and probably results 
from the mentioned above volume effect on the marginalized posterior. For the same reason the 2a region at large 
masses also becomes enhanced and the FP/HB mode again becomes strengthened. We interpret this as a reflection 
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17 (13,22) 


16.53 


16 


42% 


Non-LHC + ar + XENONioo 


122 (127,741) 


600 (608,820) 


677 (82, 1283) 


11 (9,16) 


22.21 


see text 


see text 



Table V: Best-fit points and 68% central credible regions for the CMSSM's parameters, calculated with Eq. [5} when 

constrained by different sets of experimental data. Masses are in GeV. The best-fit values of niQ are close to the 
lower edge of our prior range for niQ, which is 100 GeV; therefore, our central credible regions exclude the best-fit 

values of niQ. 



Constraint 




mi, 


B Xsj 


sin^ 0eff 


M„ 




B„ TV 






Total 




0.01 


1.32 


1.98 


4.16 


1.49 


9.76 


0.03 


0.25 


3.42 


22.21 



Table VI: Breakdown of the main contributions to for our best-fit point from a scan with a likelihood from 
non-LHC experiments and the and XENONioo limits. Note that all likelihoods are normalized to unity, 

including one for each Hj bin. 



of the fact that, despite new strong limits from the LHC, the constraining power of experimental data remains 
insufficient to eliminate prior dependence of the CMSSM. On the other hand, the effect of assuming the linear prior 
for )«i/2 and mQon the parameters tan /5 and Aq appears to be relatively weaker, as expected, since for the latter the 
prior has been assumed to be the same (linear) in both cases. 

On the other hand, the location of the best-fit point should, by definition, be entirely determined by the likelihood 
function, and therefore independent of the prior. However, because in practice it is found numerically with a Monte 
Carlo method, it is a random variable with an error and with a weak dependence on the scanning algorithm, rather 
than an exact solution. This is clear when one remembers that by choosing, for example, the log prior for the 
CMSSM's mass parameters, one in practice also chooses a log metric0 for these parameters. As a result, the 
scanning algorithm inevitably explores the low-mass region of the CMSSM's {niQ, mi/2) plane in greater detail than 
the high-mass region. This metric dependence can be overcome, or at least reduced, by tightening the algorithm's 
stopping conditions, so that it runs until it has explored the whole parameter space in sufficient detail. We also 
emphasise once again that, in Bayesian statistics the best-fit point has no significance. 

In Table |V] we present best-fit points from two scans of the CMSSM: one with non-LHC constraints only and 
one with aj and XENONioo constraints added to the likelihood. (We quote our best-fit points rounded to the 
nearest whole GeV or unit of tan /5.) In both cases the best-fit point is located in the SC/AF region of the CMSSM's 
parameter space. Our best-fit point for the non-LHC-only case is in good agreement with best-fit points reported 
previously in Ref. 0,[i3,EEE3]/ which is encouraging. Our p-value of 42% is also reasonably close to MasterCode's 
latest result for the CMSSM of 37% for this caseRef. [21]. 

On the other hand, after including the aj and XENONioo constraints in the likelihood, our best-fit point shifts 
up almost vertically to a larger value of mi/2 — 600 GeV, while the one in Ref. [21;] moves much more radically, 
to much larger values of both mass parameters (ffiQ = 450 GeV and ?Mi/2 = 780 GeV), and also giving much larger 
tan /5 = 41, although with large standard deviations reported for all the parameters. 

We investigated which constraints included in the likelihood played the most important role in determining the 
shape of the posterior, and also the location of the best-fit point. In Table |Vl] we show a breakdown of the main 
contributions to the for the best-fit point corresponding to including the ocj and XENONioo constraints in the 
likelihood. (Note that the a.j likelihood is itself a product of eight likelihoods, and therefore contributes 8 degrees 
of freedom, and inevitably has a relatively poor x^-) ft is clear that the constraint from S{g — 2)^,'^^^ plays the 
biggest role in increasing the x^ for the best-fit point. 

The upper panels of Fig.fTTlshow that the S{g — 2)^^^^ constraint requires that increases in mi/2 induced by the aj- 
limit are compensated by increases in niQ [Fig. lTTall , and also in tan fi [Fig. liibl . Our best- fit point is, however, being 
pulled in a different direction mostly by BTZ (S Xsj). The lower panels of Fig. [TTl show that increasing ?wi /2 by 
improving LHC mass limits, does not require one to increase mq [Fig. [TTcil and small tan /5 is sufficient [Fig. Iiidt , 
in order to maintain a good fit to BTZ {B Xgy). There is a clear tension between on the one hand S{g — 2)^,^^^, 
which favors lighter mass spectra in order to generate large enough SUSY contribution to the variable, and, on the 
other hand, BTZ (S —> Xgj) and other constraints which prefer a SUSY contribution to be all but suppressed]! The 
tension between the two observables is exacerbated by adding the constraint from ecj, because it pushes ?Mi/2 into 



We say metric, rather than prior, to stress that this effect is not related to our choice of Bayesian statistics. 
5 The tension between S{g — 2)^^^ and BTZ (B — )• Xs7) has already been investigated in Ref. for the pre-LHC case. 
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a region of the {mo, mi/2) plane in which it is more difficult to satisfy both constraints simultaneously. It appears 
that, at the end S{g — 2)^,^^^ is outweighed by the other constraints, most notably BTZ {B Xgj). As a result, we 
find that the best-fit point remains at small mq and tanjS, and S{g — 2)^,^^""^ is forsaken. Consequently, the best-fit 
points have a large from S{g — 2)^^^'^but not from BTZ (B Xsj); compare Table IVll 

To investigate its effect, we repeated our analysis without the S{g — 2)^^^^ experimental constraint. With a 
likelihood from non-LHC experiments, except for S{g — 2)^,^^^, and from aj, the credible regions on the (ffZO/^i/i) 
plane are similar to those in Fig. 3b though the central region of the plane is disfavored. The best-fit point's value 
of mi/2 is significantly larger than that in Fig. |3b[ This suggests that, while it favors smaller values of the SUSY 
mass parameters, S{g — 2)^^^^ does not play the dominant role in determining the credible regions. 

This is further illustrated in of Fig. Ii2ai where we show the interplay between combinations of constraints from 
BTZ [B — )• Xsj), BTZ {Bs — > fi^fi^), 3{g — 2)^^^^ and O^h^ in 2D marginalized posteriors. We can also see that 

the resulting range of BTZ (Bg — > is very close to its SM value and that O^h^ plays a rather neutral role in 

determining both the posterior and the best-fit point. This is because one can relatively easily adjust Aq to produce 
the correct value of O^h^, without much affecting the other major constraints. 

Another note is in order about ^rid p-value. As can be seen from Table El applying only non-LHC 
constraints, we find x^ = 16.53 and the p-value of 42%. Adding new constraints from the ocj- and XENONioo 
likelihood inevitably increases the value of x^/ but a shift in the p-value depends on our assumptions about the 
additional degrees of freedom. In the CMS analysis [3^] the observed events in each bin are treated as independent 
of the other bins, and described by a Poisson distribution. This is why in our treatment of the Oij limit we evaluated 
likelihood in each of the eight energy bins independently, treating them as contributing as many new degrees of 
freedom. In this case the p-value in fact increases to 61%. On the other hand, one could reasonably expect that 
signal events in different bins would be correlated, in which case the a 7- constraint could be treated as effectively 
contributing only one extra degree of freedom. In this case the resulting p-value would be 22.3%. This illustrates 
the difficulty in correctly estimating the number of degrees of freedom and the p-value. We also note that in a 
recent analysis of the MasterCode [21] the p-value diminished from 37% to 15% when the LHC and XENONioo 
constraints were included. Note, however, that in that analysis Gaussian distributions were assumed, unlike here. 
Also, their way of including LHC limits from jets+£jin the likelihood contributed only 1 additional degree of 
freedom and also forward-backward scattering variables Afl, were included in their likelihood, which suppressed 
their p-values further. All this makes it difficult to compare the resulting values of x^ arid p-value. 



In addition to a different location of the best-fit point and the p-value, the confidence intervals on the (niQ, mi/2) 
plane that are found in Ref. [21] are much larger than our confidence intervals and span the central region. On 
the other hand, they do not include the focus point region. In contrast, we find regions of the largest posterior not 
only in the SC/AF region but also in the focus point region, both in the pre-LHC case, in agreement with previous 
studies using the SuperBayeS code I14], and also after including the a.j and XENONioo limits. A physical reason 
for its existence has been given earlier. 

The apparent discrepancy between some of the results reported here and those in a recent x^ analysis by the 
MasterCode group [2r] is probably caused by a combination of several factors. Some have already been mentioned 
above. Most likely, the different ways of implementing LHC limits play an important role0 We also note that the 
implementation of BTZ {B ^ Xsj) in the version of SuperBayeS that we have used for this analysis differs from 
that of the MasterCode group whose experimental constraint is the ratio of the measured branching ratio to its 
Standard Model prediction. In contrast, in SuperBayeS the constraint from BTZ{B Xgj) is simply applied to 
the measured branching ratio. On the other hand, the choice of different statistics (Bayesian vs x^) is probably 
of secondary importance since our results for the best-fit point agree reasonably well in the non-LHC case with 
several other analyzes, both Bayesian and x^/ [i4,'i8!,[i3.Eal- We also note that the large la errors reported on the 
CMSSM parameters for the best-fit point in [25] indicate to us that probably small differences in the likelihood 
fiinction may lead to large shifts in the location of the best-fit point, in addition to numerical issues related to 
using different scanning algorithms. This may imply that the high probability regions of the CMSSM parameter 
space after including current LHC constraints are quite "unstable," or can fairly easily shift within a rather wide 
plateau, in the sense that "secondary" issues such as the precise implementation of some of the main constraints, 
or a treatment of LHC limits may lead to major differences, and seemingly contradictory results. 



^ L.R. and J. Ellis, (private communication) 
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Figure ii: Scatter plots of points from a scan with a likelihood from the ccj and XENONioo limits and from 
non-LHC experiments, colored by the values of S{g — 2)^^^"^ (top row) and BTZ (B Xs'y) (bottom row). The 
colors show the discrepancy between the CMSSM's predicted value and the central experimental value, measured 

in experimental errors. 



IV. SUMMARY 



We have performed an updated Bayesian analysis of the CMSSM. In addition to updating experimental inputs 
in indirect modes of constraining SUSY [most notably BTZ (Bg pi^jx^)], we included much improved limits from 
the CMS ocj 1.1 /fb analysis, which currently gives the strongest bounds on the CMSSM mass parameters, and from 
the XENONioo experiment. We simulated the CMS aj 1.1/fb analysis in an approximate but methodologically 
correct way by estimating the efficiency for the ocj method and constructing a likelihood function. We validated our 
method against the official CMS 95% contour. For the XENONioo limit we constructed a conservative, approximate 
likelihood function by taking into account large uncertainties related to the inputs to hadronic matrix elements and 
the local density of dark matter. 
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Figure 12: Marginalised posterior pdf for combinations of the experimental observables BTZ (S Xsj), 
3{g — 2)^,^^"^, BTZ {Bs and Oh^, from a scan with a likelihood from aj, XENONioo and non-LHC 

experiments. 



We incorporated these likelihoods into a global Bayesian fit of the CMSSM and identified marginalized posterior 
maps of the CMSSM's parameter space. These credible regions were compared with credible regions before the 
new experiments to illustrate the effects of the new constraints. The dj limit has taken a deep bite into regions 
of the CMSSM's parameter space that were previously favored, and has pushed the best-fit point to significantly 
higher values of ffii/2- Including XENONioo in the likelihood had a weak additional effect. We find that, although 
the focus point region is disfavored by the new constraints, it is still not excluded. 

Despite the disappointing null results of SUSY searches at the LHC so far, which have pushed the allowed 
ranges of CMSSM mass parameters, especially ;«i/2/ up to much larger values, we note that the best-fit point, both 
before and after including the CMS dj 1.1 /fb limit, is found close to the bottom of the marginalized la posterior 
ranges of both niQ and ?Mi/2- Despite the uncertainties in determining its location discussed above, this is certainly 
encouraging for prospects of finding a signal of SUSY in a much larger dataset already collected by both ATLAS and 
CMS. On the other hand, we note that the preference for the low-mass spectrum in the CMSSM and similar unified 
models is driven primarily by a single constraint from the 3{g — 2)^^^^^ anomaly, satisfying what has already been 
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in some tension with BTZ (B — )• Xs7), and is now becoming increasingly harder also with improving LHC limits. 
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